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Abstract. Interaction of cosmological gamma ray burst radiation with the dense 
interstellar medium of host galaxy is considered. Gas dynamical motion of inter- 
stellar medium driven by gamma ray burst is investigated in 2D approximation 
for different initial density distributions of host galaxy matter and different total 
energy of gamma ray burst. The maximum velocity of motion of interstellar medium 
is 1.8 • 10 4 km/s. Light curves of gamma ray burst afterglow are calculated for 
set of non homogeneous density, distribution gamma ray burst total energy, and 
different viewing angles. Spectra of gamma ray burst afterglow are modeled taking 
into account conversion of hard photons (soft X-ray, hard UV) to soft UV and optics 
photons. 
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1. Introduction. 

Although gamma-ray bursts (GRBs) were discovered more then thirty 
years ago (Klebesadel , 1973), their origin is still unclear. The most 
extensive data on the detection of GRBs have been obtained by the 
Compton Gamma Ray Observatory BATSE experiment (Briggs , 1995; 
Fishman , 1995; Meegan , 1992). Analysis of the GRBs detected showed 
that their apparent distribution on the sky was isotropic, but that there 
was a significant departure of their (logN - logS) curve from the N ~ 
S* -3 / 2 , low corresponding to a spatially uniform source distribution 
(Briggs , 1995; Kouveliotou , 1994). 

Observation of optical afterglows of GRB, following after identi- 
fication of GRB with a transient X-ray source by Beppo-SAX, and 
discovery of large (up to z = 4.5) redshifts in the spectra of optical 
transients had confirmed the cosmological origin of long GRB. 

The cosmological model suggests that the GRB source are located in 
distant galaxies (within ~ 10 3 Mpc). In the framework of this model, 
the observed fluxes ~ 10~ 4 erg s _1 require the release of enormous 
amounts of energy (~ 10 51 — 10 53 erg) within a fairly short time interval 
(of the order of several tens of seconds) . Such a powerful energy release 
should have a strong effect on a large volume of matter around the 
parent galaxy and should give rise to the formation of GRB counter- 
part at other wavelengths. Without specifying the mechanism for GRB 
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formation, we assume only the existence of a strong flux gamma radia- 
tion and consider the interaction of this radiation with the interstellar 
medium on large spatial and temporal scales. 

We investigate the response of an interstellar medium of standard 
chemical composition to the passage of a short-term powerful pulse 
of gamma radiation (that is, the dynamical behaviour and radiative 
cooling after heating by the gamma rays). 

The spherical symmetric model was investigated by Bisnovatyi-Kogan 
and Timokhin (1997), 2D model, permitting to study different matter 
distribution and reaction to anizotropic GRB had been studied by 
Barkov and Bisnovatyi-Kogan (2004a, 2004b) using numerical simu- 
lations by PPM method. In what follows we represent the results from 
Barkov and Bisnovatyi-Kogan (2004a, 2004b). 



2. The main equations 

We solve the system of hydrodynamic equations, describing the motion 
of matter, together with thermal processes in axially symmetric case 

% + V(p0) = 0, (1) 



d(pv r ) d(pvf, + P) 1 d(pv r v e ) 2pv1 - pvj + pv r v e ctg 9 

— 1 — r — 1 = P-fj, 



dt dr r 89 



(2) 



d(pve) d(pv r v e ) 1 d(pvj + P) 3pv r v e + pvj ctg 9 _ ^ 
dt dr r 86 r 



¥t{^ + p£ ) +v \ p \Y + £+ T)} = pH ^ pC ^ (4) 

Here p,P,e,v r , and vg are density, pressure, internal specific energy 
and two velocity components, respectively. The gamma ray pulse is 
considered as instant one having total energy T and luminosity 

L = rs (t - . (5) 

We are interested in the behaviour of the gas heated by GRB, which 
is taken as fully ionized one. 
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Consider flat spectrum of GRB 

— ^ r —E/Emax (ft\ 

dE~ E K ' 

The main GRB photons have energies larger then ionization energies 
of most electrons, so we consider energy exchange of GRB with the gas 
due to Compton and inverse Compton processes only. The function H 1 
in (4) is written us 

L He^T E max fh{,E m ax) 4/cT fc{E max ) 

7 47rr 2 m u m e c 2 

where 

fc(E max )= W(E,E max )q(E)dE, 
Jo 

fh{E max ) = - / W(E,E max )s(E)EdE. (8) 

We have used GRB spectra with E max = 0.6 MeV and 2 MeV. The 
functions s(E) and q(E), which take into account deviations from 
Thomson cross section ut due to Klein-Nishina correction {ukn\ ar e 
taken from (Beloborodov and Illarionov, 1995). They are normalized so 
that s(E) = q(E) = 1 at E < m e c 2 . In our cases f h = 0.19; f c = 0.33 
for E = 0.6 MeV, and f h = 0.065; f c = 0.16 for E = 2 MeV. 
The radiation force due to electron scattering is written as 

= ^^(U (9) 

where the function 

1 f°° 

ff(E max ) = — W(E,E max )a KN (E)dE (10) 

CT JO 

takes into account KN corrections, ff = 0.5 for E max = 0.6 MeV and 
ff = 0.32 for E max = 2 MeV. Cooling of the gas is due to different 
radiative processes (ff, fb, bb). For optically thin plasma, which is used 
for description of cooling, the function C 7 was calculated in Kirienko 
(1993), Raymond, Cox and Smith (1976) 

C^^, (11) 



volg_2004_kluw.tex; 2/02/2008; 13:45; p. 3 



4 Maxim V. Barkov and Gennady S. Bisnovatyi-Kogan 

where A(T) from Kirienko (1993) was approximated analytically with 
a precision not worse then 5% 



f °' 




T < W A K 


10" 


48.8^6.4 


10 4 < T < IO 4 25 


10" 


16.5ji-1.2 


1Q 4.25 < y < 1Q 4.5 


10" 


27.48yl.24 


1Q 4.5 < y < 1Q 5 


10" 


21.03y_o.05 


10 5 < T < IO 5 4 


io- 


13.6698y-l.413 


1Q 5.4 < T < 1Q 5.86 


10' 


22.8378y0.1515 ' 


1Q 5.86 < T < 1Q 6.19 


10" 


13.1969 ji-1.406 


1Q 6.19 < y < 10 6.83 


10" 


22. 2877^-0.075' 


1Q 6.83 < T < 1Q 7.5 


. io- 


26.6 T 0.5 


1Q 7.5 < y 



(12) 



It was shown in Barkov and Bisnovatyi-Kogan (2004a) that the heat 
conductivity may be neglected in this problem. 



3. Numerical results 

In the paper of Barkov and Bisnovatyi-Kogan (2004b) more then 10 
variants with different density distribution and GRB beaming have 
been calculated. Here we represent the main results of these calcula- 
tions. 

In the fig. 1 from Barkov and Bisnovatyi-Kogan (2004b) the evolu- 
tion of temperature distribution in the cloud with time is represented 
for the GRB exploding in the center of a spherically symmetric uniform 
cloud with a radius R = 1.5 pc, concentration nu = 10 5 cm -1 , T = 10 52 
erg, E max > 1/6 MeV. It was shown by Barkov and Bisnovatyi-Kogan 
(2004b), that for E max > 1/6 MeV, the heating at big distance from 
GRB (r > 0.05 pc) depends only on the GRB energy V and does 
not depend on E max . The temperature inversion is developed in the 
middle radiuses (r = 0.15 -r- 0.7 pc), where after heating T ~ 10 6 K, 
and the cooling is the most effective. The cooling front is propagating 
outward with superlight speed, and inward with sublight speed (phase 
velocities). The light curve for optical and ultraviolet luminosity, ob- 
served by the distant observer is represented in fig. 2 from Barkov and 
Bisnovatyi-Kogan (2004b). 

GRB heating of the cloud is most intensive in central parts and leads 
to formation of the shock wave propagating outwards. The speed of the 
shock is about 2 x 10 8 cm/s for the uniform spherical cloud, what is 
much less then the light speed. Therefore, the cloud is heated mainly by 
the light signal from GRB, and effects connected with the formation 
of central shock do not influence the integral light curve, except the 
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Figure 1. The evolution of temperature distribution in the cloud with time is rep- 
resented for the GRB burst in the center of a spherically symmetric uniform cloud 
with a radius 7? = 1.5 pc, concentration uh = 10 5 cm -1 , F — 10 52 erg, E ma x > 1/6 
MeV; curves marked by number correspond to the following time moments after 
GRB: 1) 0.76 year, 2) 1.054 year, 3) 1.103 year, 4) 1.30 year, 5) 2.40 year. 
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Figure 2. The light curve for optical and ultraviolet luminosity for same parameters, 
as in fig. 1, observed by the distant observer. Time is in years. 
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radiation in the hard X-ray band produced in close vicinity of GRB 
explosion. 

In the case, when GRB explosion takes place between two dence 
clouds, or in the cavity, produced by strong anizotropic stellar wind, 
the hydrodynamic effects may be much stronger then in the case of 
the uniform cloud, and the speed of the shock increase upto ~2x 10 9 
cm/s. 

In the fig. 3 from Barkov and Bisnovatyi-Kogan (2004a) the evo- 
lution of the velocity field is represented for the case of the density 
distribution in the cloud as n = 10 5 e — ^ ~~ 2 cos(10#) cm -3^ c i ouc [ ra _ 
dius 1.5 pc, for the energy of isotropic GRB T = 1.6 x 10 53 erg, and 
Emax > 1/6 MeV. The calculation have been performed in the region 
^ 6 ^ 7T / 10, with the condition vg = on the boundaries. The 
temperature of the heated gas depend only on the distance from GRB, 
so pressure gradient is developed inside the cavity inducing the motion 
of the matter to the axis of the cone. Collision of the flow at the cone 
axis produces a cumulative effect, leading to matter acceleration along 
the axis up to velocity ~ 2 x 10 9 cm/s Barkov and Bisnovatyi-Kogan 
(2004a). 

The accelerated matter in this case have a form of the bullet. In the 
case of the explosion in the space between two spherical clouds, the 
ejected matter should have the form of an expanding ring. 

In the case of anizotropic GRB exploded in the uniform gas cloud 
the observed light curve is different for distant observer with differ- 
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Figure 4- Light curves for collimated GRB,situated on 1 pc from center of the 
molecular cloud, time in years. Observer's situated on the line GRB - MC, case 
a); deviation from this line by the angle a = 0.1 radian, case b); by the angle 
a = 0.2 radian, case c); by the angle a — tt/2 radian, case d). 

ent angular distances from the symmetry axes. Such light curves are 
represented in fig 4. from Barkov and Bisnovatyi-Kogan (2004b). Here 
anizotropic GRB is considered with angular dependence of luminos- 
ity r(0) = 10 52 e - , = 0.1 rad, the total energy of GRB 

r to t = 2.5 x 10 49 erg, density distribution n = 10 5 e — ( r / r o) cm~ 3 , 
ro = 0.2 pc, when the explosion take place at the distance 1 pc from 
the center of the cloud. The shortest optical burst of few days with 
largest luminosity ~ 10 41 ' 5 erg/s is seen by the observer, situated at 
the symmetry axis on the continuation of the line cloud center — GRB. 

The observer on the line which is perpendicular to the symmetry 
axis is observing much longer optical afterglow (~ 1000 days), but 
with accordingly lower luminosity. All these differences are connected 
with the kinematic of the light propagation from the nonuniformly and 
nonsimultaneously heated gas cloud. 



4. Discussion 

Optical afterglow, connected with the reradiation of the GRB by the 
dense enough molecular cloud could be observed as optical transient. 
Indications that GRB explosions take place in the region of star forma- 
tion filled with dense gas clouds (Sokolov, 2001; Paczynski, 1999) make 



volg_2004_kluw.tex; 2/02/2008; 13:45; p. 7 



8 



Maxim V. Barkov and Gennady S. Bisnovatyi-Kogan 



this possibility as very probable. Observation of plato in the optical 
afterglow of GRB 030329 during a month between 64 and 94 days after 
GRB detection (Ibrahimov et. al., 2003) may be connected with such 
kind of reradiation. 

The dense molecular clouds have very low temperature, at which 
dust is formed. Estimations made by Barkov and Bisnovatyi-Kogan 
(2004b) show, that dust is evaporated by GRB pulse in the cloud near 
GRB impulse direction, or in the whole cloud by isotropic GRB with 
R ~ 1 pc, and does not influence the light curve of optical, UV and 
X-ray GRB afterglows. 
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